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ABSTRACT 

In  this  paper  we  consider  the  problem  of  using  reduced  order  dynamic  com¬ 
pensators  to  control  a  class  of  nonlinear  parabolic  distributed  parameter  sys¬ 
tems.  We  concentrate  on  a  system  with  unbounded  input  and  output  operators 
governed  by  Burgers1  equation.  We  use  a  linearized  model  to  compute  low- 
order-finite-dimensional  control  laws  by  minimizing  certain  energy  function¬ 
als.  We  then  apply  these  laws  to  the  nonlinear  model.  Standard  approaches 
to  this  problem  employ  model/controller  reduction  techniques  in  conjunction 
with  LQC1  theory.  The  approach  used  here  is  based  on  the  finite-dimensional 
Bernstein/Hyland  optimal  projection  theory  which  yields  a  fixed  finite-order 
controller. 
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1  Introduction 


In  recent  years  considerable  attention  has  been  devoted  to  the  problem 
of  using  feedback  to  control  fluid  dynamic  systems.  This  problem  is  com¬ 
plex  and  particularly  difficult  when  one  is  faced  with  phenomena  such  as 
shocks.  Moreover,  these  systems  are  governed  by  nonlinear  partial  differen¬ 
tial  equations  so  that  the  natural  state  of  the  system  is  infinite  dimensional. 
If  one  assumes  that  “full  state  feedback”  is  necessary  to  design  practical 
controllers,  then  one  would  conclude  that  feedback  control  of  fluid  dynamic 
system  is  “not  practical”.  However,  it  is  well  known  that  even  in  finite  di¬ 
mensional  control  systems  one  rarely  has  the  ability  to  accurately  sense  all 
states,  so  that  some  form  of  dynamic  compensation  must  be  used. 

This  idea  clearly  extends  to  infinite  dimensional  problems  and  there  is 
a  growing  literature  on  observers/compensators  for  distributed  parameter 
systems.  In  this  paper  we  consider  a  boundary  control  problem  governed  by 
Burgers’  equation.  We  selected  this  problem  because  Burgers’  equation  is 
an  infinite  dimensional  model  that  captures  some  phenomena  (e.g.,  shocks) 
often  observed  in  fluid  flows  and  because  it  is  simple  enough  to  provide 
real  insight  into  the  problem.  The  goal  is  to  show  that  it  is  possible  to 
use  modern  control  theory  to  produce  practical  finite  dimensional  dynamic 
compensators  for  boundary  control  of  nonlinear  partial  differential  equations 
of  the  type  that  occur  naturally  in  fluid  dynamics. 

We  shall  present  a  short  summary  of  one  approach  (the  optimal  projec¬ 
tion  method  due  to  Bernstein  and  Hyland)  and  show  how  this  approach  can 
be  used  in  conjunction  with  standard  numerical  schemes  to  produce  a  real¬ 
izable  low  order  controller.  The  optimal  projection  method  is  one  of  many 
approaches  to  this  problem.  However,  we  shall  concentrate  on  this  method 
because  a  very  nice  theory  has  already  been  developed  (for  bounded  input 
and  output  operators)  and  we  are  more  interested  in  illustrating  (to  non¬ 
experts)  that  recent  results  in  distributed  parameter  control  theory  can  be 
used  to  design  practical  feedback  laws,  than  in  discussing  the  “best”  ap¬ 
proach  to  the  problem.  It  will  be  clear  from  our  presentation  that  we  are 
writing  for  those  that  are  not  necessarely  “control  experts” .  The  extension 
of  the  general  theoretical  results  to  unbounded  input  and  output  operators 
will  appear  in  a  forthcoming  paper.  However,  for  the  compensators  pre¬ 
sented  here,  we  do  not  need  the  most  general  theory  since  we  use  the  finite 
dimensional  version  of  the  optimal  projection  method. 

As  noted  above  it  is  almost  impossible  to  observe  the  whole  state.  Con- 
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trols  and  sensors  are  limited  to  a  few  points  or  segments  of  the  boundary,  so 
it  is  necessary  to  construct  an  appropriate  observer  (estimator)  of  the  state 
and  design  a  feedback  control  law  (called  a  compensator)  based  on  the  in¬ 
formation  available  from  the  observed  (estimated)  state  variable.  Boundary 
control  and  observation  often  leads  to  unbounded  input  and  output  opera¬ 
tors.  Stabilization  by  dynamic  feedback  or  compensation  has  been  consid¬ 
ered  by  Curtain  [5],  Fujii  [8],  and  Nambu  [12]  for  classes  of  parabolic  as  well 
as  hyperbolic  systems,  including  control  and  observation  at  the  boundary. 
All  of  these  approaches  produce  stabilization  schemes  that  cither  have  the 
same  finite  order  as  that  of  a  high-order  approximate  model,  or  alternatively, 
open-loop  model  reduction  or  closed-loop  control  reduction  techniques  are 
applied  to  achieve  a  lower-order  compensator.  An  advance  was  made  by 
Schumacher  [15],  when  he  gave  a  theory  for  designing  finite-dimensional 
compensators  for  a  large  class  of  systems,  including  parabolic  and  delay  sys¬ 
tems.  However,  in  his  theory  it  was  assumed  that  the  control  and  observation 
operators  are  bounded.  Curtain  [4]  presented  an  alternative  compensator 
design  which  applied  to  the  same  class  of  systems,  except  that  unbounded 
inputs  and  outputs  were  allowed.  In  [14],  Pritchard  and  Salamon  established 
a  framework  based  on  semigroup  theory  for  treating  the  linear  quadratic  reg¬ 
ulator  problem  for  infinite-dimensional  systems  with  unbounded  input  and 
output  operators.  Their  approach  is  based  on  a  weak  formulation  of  the  Ric- 
cati  equations  which  characterize  the  optimal  feedback  law  in  an  appropriate 
dual  space. 

Here  we  consider  the  problem  of  designing  a  fixed-finite-dimensional  com¬ 
pensator  for  a  class  of  distributed  system  governed  by  Burgers’  equation, 
where  the  control  and  the  observation  are  implemented  at  the  boundary  of 
the  domain.  The  possibility  of  applying  this  approach  to  distributed  pa¬ 
rameter  systems  was  first  suggested  by  Johnson  in  [9]  and  Pearson  [13]. 
The  idea  of  fixing  the  order  of  the  finite-dimensional  compensator,  while 
retaining  the  distributed  parameter  model  was  expanded  and  developed  by 
Bernstein  and  Hyland  in  [1]  and  [2].  The  method  extends  the  full  order 
LQG  case  to  an  “optimal  fixed-finite-order  compensator”  characterized  by 
four  equations;  two  modified  Riccati  equations  and  two  modified  Lyapunov 
equations,  coupled  by  an  oblique  projection  whose  rank  is  precisely  equal  to 
the  order  of  the  compensator.  Bernstein  and  Hyland  assumed  that  the  con¬ 
trol  and  observation  operators  were  bounded  and  hence  boundary  control 
and  observations  were  not  covered  by  their  theory. 

We  will  present  a  Bernstein/Hyland  type  fixed-finite-dimensional  com¬ 
pensator  design,  which  does  extend  to  unbounded  input/output  problems. 
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In  Section  2  we  discuss  the  existence  of  a  finite-dimensional  compensator  for 
parabolic  distributed  parameter  systems  with  unbounded  control  and  obser¬ 
vation.  In  Section  3  we  summarize  the  infinite-dimensional  optimal  projec¬ 
tion  theory  from  [1],  and  derive  the  corresponding  equations  and  feedback 
gains  which  characterize  the  fixed-finite-order  compensator.  In  Section  4  we 
present  an  example,  construct  the  approximation  schemes,  and  discuss  the 
computational  algorithm  used  for  the  optimal  projection  design  synthesis. 
Finally,  Section  5  contains  numerical  results  and  Section  6  is  devoted  to  a 
few  closing  remarks. 

2  A  Theoretical  Existence  Result 

We  consider  the  following  abstract  Cauchy  problem 

z(t)  =  Az(t)  +  Bu(t),  z(0)  =  z„  G  H  (1) 

y(t)  =  Cz(t)  t>  0  (2) 

where  H  is  a  Hilbert  space,  u(.)  G  £2( 0  y(.)  G  L2(0  T;S^),  and 

A  is  the  infinitesimal  generator  of  analytic  semigroup  S{t)  on  H,  generally 
unstable,  with  exponential  growth  rate 

W°  =  fe*"1  l°gll5(*)ll£(ff)  >  0  (3) 

so  that 

||S(t)||£(ff)  <  Me for  ail  e  >  0,  t  >  0  (4) 

for  some  constant  M  =  M{wo,e)  >  1.  Throughout  the  remainder  of  this 
paper  we  let  A  denote  the  translation  A  =  —A  +  wl ,  where  w  is  fixed 

and  w  >  wo,  so  that  A  has  well-defined  fractional  powers  (A),J  on  H  and 

—A  is  the  generator  of  a  strongly  continuous  analytic  semigroup  S(t)  on  H 
satisfying 

\\S(t)\\c{H)  <  Me"*',  t>0.  (5) 

In  order  to  allow  for  unbounded  operators  B  and  C,  we  assume  that  B 
G  £(5?m,  V )  and  C  G  C(W,  3^),  where  W  and  V  are  also  Hilbert  spaces  such 
that 

V(A)  CW'-tH’-iV  (6) 

with  continuous  dense  injections.  More  precisely,  we  assume  that  B *  is 
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|A*]7-bounded,  or  equivalently, 

[A]_7B  €  C(%tm,H)  for  0  <  7  <  1.  (7) 

Similarly,  for  the  operator  C  we  assume  that 

C[i]'7  g  C(H,&)  for  0  <  7  <  1.  (8) 

It  is  helpful  to  interpret  (1-2)  in  mild  form.  In  particular,  the  solution 
2(1!)  is  given  by 

z(t)  =  S(t)z0  +  f  S(t  —  s)Bu(s)ds.  0  <t<T  (9) 

Jo 

and  the  output  by 

y(t)  -  CS(t)zo  +C  f  S(t  -  s)Bu(s)ds.  (10) 

Jo 

We  assume  that  S(t)  is  also  an  analytic  semigroup  on  W  and  that  the 
following  hypotheses  are  satisfied: 

(H-l)  There  exists  a  constant  b(T)  >  0  such  that  for  every  T  >  0, 
Jq  S(T  —  s)Bu(s)ds  €  W  and 

S(T-s)Bu{s)ds^  <  6(r)||«(.)||t»(o  t;»-)  (ID 

for  every  u(-)  6  L2( 0  T;  3?m). 

(H-2)  There  exists  a  constant  c(T)  >  0  such  that  for  every  T  >  0, 

£  ||CS(t)x|||3(0  r.^dt  <  c(T)||x||v  (12) 

for  every  x  €  W . 

We  now  give  sufficient  conditions  which  imply  that  the  system  (1-2)  can 
be  stabilized  by  a  finite-dimensional  compensator  of  the  form 


w(t)  =  Acw(t)  -  Bcy(t) 
u(t)  =  Ccw(t ) 


w(0)  =  w0 


(13) 

(14) 
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where  Ac  €  3tN‘xNc,  Bc  G  %tNcXm,  and  Cc  €  $texNc  are  suitably  chosen 
matrices.  We  need  the  following  well-posedness  result  for  the  connected 
system  (1-2)  and  (13-14).  This  result  and  proof  may  be  found  in  [4]. 

Proposition  2.1  Let  (H-l)-(H-2)  be  satisfied,  then  for  all  zQ  G  W ,  wa  G 
there  exists  a  unique  solution  pair  z(t)  and  w(t )  of  (1-2)  and  (13-14)- 
This  means  that  z(t)  is  continuous  in  H  and  absolutely  continuous  in  V, 
that  (1)  is  satisfied  for  almost  every  t  >  0  where  u(t)  is  given  by  (14),  and 
that  w(t)  G  %lNc  is  continuously  differentiable  and  satisfies  (13)  where  y(t) 
is  given  by  (2). 

In  addition  to  hypotheses  (H-l)  and  (H-2),  we  assume: 

(H-3)  Stabilizability  Condition  (S.C.) 

There  exists  an  operator  F  G  C(H,  3Jm)  such  that  Af  =  A  +  BF 
generates  an  analytic  semigroup  SF(t)  =  e^A+BF^‘  and  SF(t)  is  expo¬ 
nentially  stable  on  H,  i.e., 

||SF(*)!!£(ff)  <  MFe~WFt,  for  wF  >  0.  (15) 

(H-4)  Detectability  Condition  (D.C.) 

There  exists  an  operator  G  G  £(3^,  H)  such  that  Ag  —  A  +  GC 
generates  an  analytic  semigroup  Sc{t)  =  e(-A+GC^t  and  So(t)  is  expo¬ 
nentially  stable  on  H,  i.e., 

ll5G(t)ll£(//)  <  Mce_u,Gi,  for  wG  >  0.  (16) 

(H-5)  In  addition  to  (H-3)  and  (H-4)  there  exists  a  finite-dimensional 
subspace  R  C  W,  with  dimR  <  Nr  such  that 

(i)  S/r(t)R  c  R,  for  all  t  >  0, 

(ii)  Range  GCK, 

(iii)  R  C  V(Af). 

Moreover,  there  exist  linear  maps  i:  3?‘Vc  — >  R,  n  :  H  — »  3?  v<r  such  that 
7 rt  = /atc,  inx  =  x  for  x  G  R.  (17) 

Note  that  (H-5)  implies  that  xAFi  is  a  well  defined  linear  map  on  f?'Nc.  We 
will  show  that  the  system 

w(t)  =  7r  (Af  +  GC)  i  w(t)  —  xGy{t),  w( 0)  =  wa  (18) 
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u(t)  =  Fiw(t)  (19) 

defines  a  stabilizing  compensator  for  the  Cauchy  problem  (1-2).  The  follow¬ 
ing  result  is  a  slight  extension  of  Theorem  2.5  in  [7]  for  unbounded  inputs 
and  outputs. 

Theorem  2.2  If  (H-l)-(H-5)  are  satisfied,  then  the  closed-loop  system  de¬ 
fined  by  (1-2)  and  (18-19)  is  exponential  stable. 

Proof:  Note  that  without  loss  of  generality  we  can  assume  that  dim  N  =  Nc. 
By  Proposition  2.1  it  follows  that  the  closed-loop  system  is  a  well-posed 
Cauchy  problem.  Let  z0  €  VP,  w0  £  and  z(t),  w(t)  be  defined  by  (1-2) 
and  (18-19),  respectively.  Since  z(t)  £  VP,  if  x(t)  is  defined  by 

x(i)  =  iw(t)  —  z(t)  t  >  0, 

then  x(t)  belongs  to  VP  and  it  is  straightforward  to  show  that 

w{t)  —  it  AFiw(t) +  ttGC  x(t).  (20) 


Therefore, 


=  inSp(t)iw0  +  f  in Sp(t  —  s)in GC x(s)  ds  —  z(t) 
JO 

—  SF(t)iw0+  I  Sf(t  —  s)  GC  x(s)  ds  —  z(t) 

Jo 

—  S(t)iw0+  f  S(t  —  s)[B  F  iw(s)  +  GC  x(s)]ds 

Jo 

—  S(t)zo+  f  S(t  —  s)  B  u(s)  ds 
Jo 

=  S(t)x( 0)+  /  S(t-  s)GCx(s)ds, 

Jo 


which  implies  that  x{t.)  =  Sg'(£)z(0).  The  stability  of  x(t ),  w(t)  and  z(t.) 
follows. 
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3  Optimal  Projection  Theory 

Consider  the  steady-state  fixed-order  dynamic  compensator  problem,  de¬ 
fined  by  the  infinite-dimensional  control  system 

z(t)  =  Az(t)  +  Bu(t)  +  HiTj(t)  (21) 

with  measurements 

y(t)  =  Cz(t)  +  H2y(t).  (22) 

The  objective  is  to  design  a  finite-dimensional  fixed-order  dynamic  compen¬ 
sator 


zc(t)  =  Aczc(t)  +  Bcy(t)  .(23) 

u(t)  =  Cczc(t )  (24) 

which  minimizes  the  steady-state  performance  criterion 


J(AC,  Bc,Cc)  ==  lim  7  /  E  (Riz(s) ,  z(s))  -f  u(s)T  R2u{s) 
<->oo  t  Jo 


ds  (25) 


where  the  operators  A,  B  and  C  satisfy  all  the  assumptions  given  in  the 
previous  section  and  E[  ]  is  the  expectation.  In  addition,  asstime  that  the 
state  and  measurements  are  corrupted  by  a  white  noise  y(t)  in  the  Hilbert 
space  H,  with  zero-mean  Gaussian,  Hi  €  C(H,  H ),  H2  £  C(H ,  5R*),  R\  £ 
C(H)  is  self-adjoint  and  nonnegative  definite,  and  that  R2  is  an  m  x  m 
symmetric  positive-definite  matrix.  We  assume  that  the  disturbance  and 
measurements  are  independent,  i.e.,  Hi  H2  =  0,  Vj  =  Hi  H\  6  C(H)  is 
nonnegative  definite  and  of  trace  class,  and  that  V2  —  H2H2  £  is 

positive  definite.  Also,  it  is  assumed  that  the  initial  state  z(0)  =  z0  is 
Gaussian  and  independent  of  y(-).  The  compensator  will  be  assumed  to  be 
of  fixed,  finite  order  Nc  (i.e.,  zc(t )  £  RNc)  and  the  optimization  is  performed 
over  Ac  £  3?*'**',  Bc  £  and  Cc  £ 

If  one  introduces  the  augmented  state  space  ~H  —  H  x  RN-,  then  the 
closed-loop  system  becomes  a  linear  system  on  H.  Consequently,  define  the 
closed-loop  operator  A  :  T>(A)  C  H.  — >  H  on  the  dense  domain  X>(*4)  = 
V(A)  x  &N'  by 


A  BCc 
BcC  Ac 


‘  A 

0  ' 

+ 

0 

BCc 

0 

°. 

BcC 

Ac 
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Since  the  operator 


A  0 
0  0 


generates  an  analytic  semigroup 


eAt  0 

0  7„c 


t  >  0. 


then  conditions  (7)-(8)  imply  that  A  is  also  closed  and  generates  an  analytic 
semigroup  e*4*  on  (see  [10]).  To  guarantee  that  J  is  finite  and  indepen¬ 
dent  of  initial  conditions  we  restrict  our  attention  to  the  set  of  admissible 
compensators  defined  by 

S  =  {(.i4c,  Bc,  Cc)  :  eAt  is  exponentially  stable}.  (26) 


Tf  ( Ac ,  Dc,  Cc)  €  S ,  then  there  exist  a  >  1  and  /3  >  0  such  that 

\\eA,\\  <  ae~0t  t  >  u.  (27) 

Moreover,  we  know  from  Theorem  2.2  above  that  S  is  non-empty.  We  now 
state  some  results  found  in  [1]  and  [2], 

Lemma  3.1  If  Q  and  P  €  C(H)  have  finite  rank  and  are  nonnegative 
definite,  thenQP  is  nonnegative  semi-simple.  Furthermore,  if  rank  ( QP )  = 
Nc,  then  there  exist  G  and  T  €  £(//,5?'c)  and  a  positive  semi- simple  matrix 
M  such  that 


QP  =  G*Mr  (28) 

rG*  =  ISc.  (29) 


Proof:  Bernstein  and  Hyland  give  a  complete  proof  of  this  result,  in  [1].  Here 
we  outline  their  proof  in  order  to  illustrate  the  form  of  the  factorization  of 
QP  and  to  provide  a  description  of  the  operators  G  and  T.  Since  Q  and 
P  have  finite  rank,  there  exists  a  Unite  dimensional  subspace  Z  C  H  such 
that  QZ  C  Z,  QZ1  =  0,  PZ  C  Z  and  PZ 1  =  0.  Hence  there  exists  an 
orthonormal  basis  for  H  and  in  this  basis  Q  and  P  have  the  infinite  matrix 
representations 


Q  i  o 
0  0 


p  = 


p\  o 
0  0 
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where  C^i,  Pi  €  9?rXr  and  r  =  dim  Z.  Consequently,  there  exists  an  invert¬ 
ible  $  G  3?rxr  such  that  A  =  $-1QiPi$  is  nonnegative  and  diagonal  and 
QP  is  nonnegative  and  semi-simple.  If  rank  (QP)  =  Nc ,  then  it  is  clear  that 
$  can  be  chosen  so  that 


where  A  G  3?‘VcX  Vc  is  positive  and  diagonal.  Hence, 


QP  = 


$  0 
0  Ia 


Ikc 

0 

0 


A  [  [  INc  0  ]  O' 


0  ' 

0  Ioc 


and  if  we  define  G ,  M  and  T  by 


G 

r 

M 


[K  0] 

[[*"  °] 


<&T  0  ' 

0  loo 

0 

0  loo 


5~*  AS. 


for  any  invertible  S  G  ?RNc*Nc,  then  G,  T  and  M  provide  the  desired  factor¬ 
ization  and  this  completes  the  proof. 


Throughout  the  paper  we  will  refer  to  G,  T  and  M  satisfying  the  above 
lemma  as  a  (G  --  M  —  T)  -  factorization  of  QP.  For  convenience  we  define 
E  =  BR^]B*  and  E  =  C'V^lC  and  let  /,vc  and  In  denote  respectively 
the  Nc  x  Nc  identity  matrix  and  the  identity  operator  on  H ,  respectively. 
We  state  Bernstein’s  and  Hyland’s  main  theorem  which  provides  a  set  of 
necessary  conditions  that  characterize  the  optimal  steady-state  fixed  order 
dynamic  compensator  for  bounded  input  and  output  operators  (see  [1]). 

Theorem  3.2  Let  B  and  C  be  bounded  operators  and  let  Nc  be  given  and 
suppose  that  there  exists  a  controllable  and  observable  dynamic  compen¬ 
sator  ( Ac ,  Bc,Cc)  G  S  of  order  Nc  which  minimizes  J  given  by  (25),  then 
there  exist  nonnegative  definite  operators  Q,  P,  Q,  and  P  on  H  such  that 
Ar,  Be,  and  Cc  are  given  by 

Ac  =  ?(A  -  QE  -  EP)G*  (30) 
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Bc  =  r  QCmVf'  (31) 

Cc  =  -R^lBmPTm  (32) 

for  some  (G  -  M  —  T)  -  factorization  of  QP  and  such  that,  with  r  =  G"T  € 
C(H),  the  following  conditions  are  satisfied: 

Q  :  V(A’)^V(A)  P  :  V{A)  ->V(Am) 

Q  :  H  — ►  T>(A)  P  :  H-^V(A’) 

rank  (Q)  =  rank  (P)  =  rank  (QP) 

and 

0  =  (A  -  tQZ)Q  +  Q(A  -  rQE)*  +  Vi  +  rQEQr*  (33) 

0  =  (A-YiPryP  +  P(A-1:Pt)  + Ri +t"PEPt  (34) 

0  =  [(A  -  DP)Q  +  Q(A  -  DP)'  +  QDQ\  r*  (35) 

0  =  [(A-QD)'P  +  P(A  -QD)  +  PDP]r.  (36) 

Note  that  these  necessary  conditions  consist  of  a  system  of  four  operator 
equations,  including  a  pair  of  modified  Riccati  equations  and  a  pair  of  mod¬ 
ified  Lyapunov  equations  which  are  coupled  by  the  operator  r  €  C(H).  The 
operator  r  is  idempotent,  since  r2  =  tt  —  G’TG'T  =  G'I„cT  =  GT  =  r. 
In  general  t  is  an  oblique  projection  and  may  not  be  orthogonal  since  there 
is  no  requirement  that  t  be  self-adjoint.  Moreover,  we  note  that  in  view  of 
Lemma  3.1,  Theorem  3.2  applies  to  ( SAcS~l,SBc ,  CcS~x )  for  any  invertible 
S  €  &NcX*c,  since  the  (G  —  M  —  Tj-factorization  of  QP ,  used  to  determine 
Ac,  Bc  and  Cc,  is  not  unique.  However,  the  operator  r  remains  invariant 
over  the  class  of  factorizations.  An  easy  computation  yields  the  following 
identities: 

Q  =  tQ  and  P  =  Pt.  (37) 

It  is  helpful  to  have  an  alternative  form  of  the  optimal  projection  equa¬ 
tions  to  actually  compute  the  optimal  fixed-order  compensator  of  the  ap¬ 
proximating  finite-dimensional  plant.  The  following  result  for  bounded  input 
bounded  output  operators  may  be  found  in  [1], 

Proposition  3.1  If  B  andC  are  bounded,  then  the  optimal  projection  equa¬ 
tions  (33)-(36)  are  equivalent,  respectively,  to 

0  =  AQ  +  QA'  +  V i  -  QDQ  +  rxQE(?r*  (38) 
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0  =  AmP  +  PA  +  Ri  -  PDP +  tIPY.Ptl  (39) 

0  =  APQ  +  QAP  +  QDQ  -  tlQT,Qt1  (40) 

0  =  A*P  4-  PAg  +  PDP  -  t±PDPtl  (41) 

where 

tl  —  Iff  -  r,  Ap  =  A  -  DP  and  Aq  =  A  -  QD.  (42) 


This  form  of  the  optimal  projection  equations  shows  that  there  is  a 
connection  between  Theorem  3.2  and  the  standard  LQG  result  when  dim 
H  =  N  <  oo.  In  this  case,  we  note  that  the  (G  —  M  —  T)-factorization  of 
QP  when  Nc  =  N  is  given  by  G  =  T  =  Is  and  M  =  QP.  Since  r  =  7,v 
and  rx  =  0,  it  follows  that  (38)-(39)  reduce  to  the  standaid  observer  and 
regulator  Riccati  equations. 

To  obtain  a  geometric  interpretation  of  the  optimal  projection  we  intro¬ 
duce  the  “quasi-full-state”  estimate 

z(t)  =  G'zc(t)  6  H,  (43) 

so  that  r z(t)  =  z(t)  and  zc(t)  =  Fz(t).  Hence,  the  closed-loop  system  can 
be  written  as 

z{t)  =  Az(t)  +  BCcTz(t)  (44) 

z(t)  =  t(A  + BCc  -  BcC)rz(t)  +  TBcCz{t)  (45) 

where 

Bc  =  QC’V2~l  and  Cc  =  -R^lB‘P.  (4G) 

This  shows  that  the  geometric  structure  of  the  quasi-full-order  compen¬ 
sator  is  dictated  by  the  projection  r.  Sensor  inputs  tBcCz  are  annihilated 
unless  they  are  contained  in  TZ(r’)  =  A^(r)1,  while  tz  employed  in  the  con¬ 
trol  input  is  contained  in  7£(r).  Consequently,  7 Z(r)  and  TZ(t')  are  the 
control  and  observation  subspaces  of  the  compensator,  respectively.  In  order 
to  modify  the  previous  results  so  that  they  will  apply  directly  to  unbounded 
B  and  C  operators,  care  must  be  exercised  to  precisely  define  the  weak 
forms  of  (33)-(36)  and  (38)-(41).  We  shall  not  consider  this  problem  in  this 
short  note.  However,  we  shall  use  these  systems  to  guide  the  approximations 
below. 
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4  Finite  Dimensional  Approximation 

In  general,  the  optimal  projection  equations  (38)-(41)  are  infinite  dimen¬ 
sional  operator  equations.  To  actually  use  these  equations  to  compute  the 
optimal  fixed-finite-order  compensator,  a  finite  dimensional  approximation 
is  needed  (see  [2]  for  details). 

Let  Hn  for  N  —  1,2,  ••*,  be  a  sequence  of  finite  dimensional  linear 
subspaces  of  H  and  let  VN  :  H  — >  HN  be  the  canonical  orthogonal 
projections.  Let  AN  €  C(HN),  BN  €  £(ftm,  HN),  CN  €  C(HN ,  9?*), 
R f  €  C(HN)  and  VXN  €  C(HN)  be  given  and  consider  the  approximating 
system 

zN(t)  =  ANzN(t)  +  BNuN(t)  +  H?r)N(t)  (47) 

yN(t)  =  CNzN(t)  +  H?VN(t).  (48) 

The  goal  is  to  design  a  sequence  of  finite-dimensional  dynamic  compensators 
of  fixed  order  Nc  of  the  form 

i"(t)  =  A?z?(t)  +  B?yN(t)  (49) 

uN(t)  =  CcNzN(t),  (50) 

which  minimizes  the  performance  criterion 

JN(A?,B?,C?)d^  Hm  jjf  B^/W,*"^))  +  u(s)TK2u(s)]  ds. 

0  (51) 

Now,  for  each  N  =  1,2,---,  let  j  be  a  basis  for  HN .  Also,  for  any 

linear  operator  FN  with  domain  and  range  in  HN ,  unless  otherwise  noted, 
we  use  the  same  symbol  FN  for  its  matrix  representation  with  respect  to 
the  basis  chosen.  Let  'I'^  denote  the  kN -square  Gram  matrix  corresponding 

to  the  basis  (e.g.,  VN  =  [(<£,*  ,  )ffW]  )•  Note  that 

(Aa’)*  =  (*N)-l(AN)T*N  (Bn)'  =  (BN)T<1>"  (52) 

(CNy  =  ('PA,)"1(Cyv)T  {EN)  =  BNR^](BN)Jy iN  (53) 

(r£)“  =  (tfWi)1**  (DN)  =  (*n)-'(Cn)tV7-xCn  (54) 
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and  if  we  define  the  kN  x  kN  nonnegative  definite  matrices 


Q6V 

def 

^(^r1 

pJV  def  yNpN 

def 

pN  def  q/N  pN 

V'o" 

def 

^($iV)-l 

R$  d4f 

def 

Dq  d^f  (Cn)tV2-'Cn 

then  the  matrix  equivalence  of  the  operator  equation  (38)-(41)  become 


0  = 

ANQ$  +  Qo(As)J  +  V0N  - 

+tlQq^0  Qo  (rx  )T 

(55) 

0  = 

WfPf  +  PfA"  +  <  -  P6v^Pjv 

+(t*)tP0;v£0vP0at'v 

(56) 

0  = 

<<?£  +  Q6v(<)t  +  Q6v^Qon 

~t?QoZXQo(t?)T 

(57) 

0  = 

(A£)tP0n  +  Pf<  +  po^o  Po 

(-N\t  pNyN  pN  N 
~\Tx  )  * o  "o  Po  Tx  • 

(58) 

The  approximating 
Nc  is  then  given  by 

optimal  dynamic  compensator  (A^,B^ 

,  C?)  of  order 

=  Tq  {AN  -  QqDq  -  Dq  Pq  )(Gq  )t 

(59) 

B* 

=  poQo(^N)rV2~l 

(60) 

C” 

=  -pj1(P,v)1>ov(ro  )T 

(61) 

where  Tq  ,  Gq  €  J?jVcXfc,V  and  Mq  £  provide  a  (Gq 

i 

f-4 

1 

1 

factorization  of  Qq  Pq  . 


We  turn  now  to  an  example.  Consider  Burgers’  equation,  with  Neumann 
boundary  control  given  by 


d  ,  > 

atz(Ux) 

z( 0,  x ) 


0_ 

dx 


z(f,0) 


Cf  d 

0<I<1,  t>  0  (62) 

z0{x)  (63) 


~«i(<), 


(64) 


—  z(t,  1)  =  u2(t), 
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and  observations 


yi(t)  =  z(t,0) 

(65) 

y*  (t)  = 

(66) 

where  e  =  >  0  and  Re  is  the  Reynolds  number.  Initially, 

the  linearized  Neumann  boundary  control  problem 

we  consider 

= 

d2 

e-^z(t,i),  0  <  x  <  1,  t>0 

(67) 

2(0,1)  = 

zo(x) 

(68) 

^-z(t,l)  =  u2(t). 

(69) 

We  will  apply  the  linearized  feedback  control  laws  constructed  from  this 
model  to  the  nonlinear  Burgers’  equation.  System  (67)-(69)  can  be  placed 
into  the  standard  state  space  framework  by  defining  the  operator  A.  on 
H  =  L2( 0,1)  by 

A(tf>  =  e<f>"  (70) 

for  all  4>  €  V(A()  =  {<f>  €  H2(0, 1)  :  =  <^'(1)  =  0}.  Define  W  =  F*  = 

H$  (0, 1)  =  V{a!  )  and  let  B  :  fi2-+V  be  defined  by  B  =  AAf  where  A  = 
— Af+wl  and  we  assume  that  w  is  not  an  eigenvalue  of  A(  with  homogeneous 
Neumann  boundary  conditions,  so  that  A  is  boundedly  invertible  on  £2(0, 1). 
The  Neumann  map  Af  is  defined  by  the  boundary  system  given  in  [11,  pages 
53-56].  Let  C  :  W  — *  5R2  defined  by 


c<t>  = 


m 

mu 


(71) 


The  boundary  control  problem  (67)-(69)  can  be  represented  by  a  differential 
equation 


=  Atz{t)  +  Bu,  z(0)  =  z0  (72) 

y(t)  =  Cz(t).  (73) 


It  is  well  known  that  A(  generates  an  analytic  semigroup  S(t)  on  H .  More¬ 
over,  the  spectrum  o(A^)  of  A(  consists  of  all  eigenvalues  An,  n  =  0, 1,2,  •  •  ■ 
given  by  An  =  — en27r2  and  for  each  eigenvalue  An  the  corresponding  eigen- 
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function  <f>„  is  given  by 

<f>o(x)  =  1  <f>n{x)  =  \/2cos(n7rx),  0  <  i  <  1.  (74) 

One  can  easily  verify  conditions  (7)-(8),  by  taking  7  =  5  (see  [11]).  It  is 
straightforward  to  show  that  (H-l)-(H-5)  are  valid. 


5  Numerical  Results 


Now,  we  formulate  a  specific  approximation  scheme  for  the  boundary  control 
problem  (74).  For  each  N  =  2, 3,  •  •  •  let  divide  the  unit  interval  [0,1]  into  N 
equal  subinterval  [x<,  x^+i],  x,  =  i  =  1, 2,  •  •  • ,  N  +  1.  Let  HN  =  Span 

where  h^(-)  are  the  standard  hat  functions  defining  continuous 

piecewise  linear  splines  (see  [3]).  Note  that  kN  =  dim  HN  =  N  +  2  and  let 
the  approximate  solution  zN (t,x)  of  z(t,  x)  for  equation  (72)-(73)  be  given 
by 

N+l 

zN(t,x)=  Y,  z?(t)f li'(x)  (75) 

i=0 


for  some  z/v(t)  €  St,  i  =  0, 1,  •  •  • ,  N  +  1.  Standard  finite  element  approxima¬ 
tions  yield  the  ODE  system 


jt^N(t)  =  A*'zs(t)  +  BNu(t),  zN(0)  =  zg  (76) 

yN(t)  =  CNzN(t)  (77) 


where  the  matrices  A? ,  BN ,  CN  can  be  easely  computed  by  using  the 
Ritz-Galerkin  approximation. 

For  our  numerical  example,  we  set  a  =  L,  the  initial  condition  z0(x)  = 
sin(7rx),  n  =  v\  =  1  and  r2  =  v2  =  10-3.  Also,  R\  =  r\IH ,  R2  =  r2Im, 
V\  —  v\ In,  and  V2  =  v2I(.  Therefore,  it  follows  from  Section  3  that  R°  = 
and  Vj*  =  ui^^)-1  where  is  the  Gram  matrix.  In  this  numerical 
example  we  will  compare  the  approximating  optimal  LQG  (i.e.,  Nc  =  N  +  2) 
with  the  dynamic  compensators  of  various  order  Nc.  The  optimal  projection 
equations  (55)-(58)  were  solved  using  the  homotopic  continuation  algorithm 
described  in  [16].  The  approximating  controllers  defined  by  the  linear  fixed- 
order  compensator  ( B ^  and  )  were  applied  to  Burgers’  equation  (62)- 
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We  note  that  N  =  32  produces  converged  optimal  LQG  designs.  Hence, 
the  reduced  order  compensators  were  tested  on  both  the  linear  and  nonlinear 
problem  using  the  N  =  32  order  finite  element  model. 

In  the  full  order  case  N  =  32  and  Nc  =  34,  the  converged  feedback  and 
observer  functional  gains  are  given  in  Figures  1  and  2,  respectively.  Since 
we  are  controlling  the  flux  at  each  end  point  x  =  0  and  x  =  1,  we  have 
two  feedback  functional  gains,  the  one  plotted  with  solid  line  is  the  flux 
control  gain  at  the  origin  and  the  one  plotted  with  dashed  line  is  the  flux 
control  at  the  end  point  x  =  1.  Similarly,  since  we  are  sensing  the  flow 
at  the  origin  and  the  end  point,  we  have  two  observer  gains  (solid  line  for 
observer  gain  at  the  origin  and  dashed  line  for  the  observer  gain  at  x  =  1). 
Next,  we  applied  the  full  order  controller  to  Burgers’  equation  resulting  in 
the  nonlinear  closed-loop  trajectory  given  in  Figure  3. 

In  the  fixed-order  case,  we  considered  the  accuracy  of  the  impulse  and 
step  responses  of  the  various  reduced  order  compensator  designs  compared 
to  the  corresponding  responses  of  the  full  order  LQG  design.  Figure  4  il¬ 
lustrates  the  linear  closed-loop  impulse  response  for  the  full-order  LQG  and 
reduced  order  compensator  (of  order  Nc  —  16)  designs.  The  impulse  re¬ 
sponse  of  the  linear  closed-loop  system  for  the  16th-order  compensator  is 
in  perfect  agreement  with  the  LQG  response.  Note  that  in  Figure  4  we  see 
only  one  plot  for  both  designs  because  both  plots  are  essentially  the  same. 
Similar  trends  are  seen  (Figure  5)  in  the  comparisons  of  the  step  responses 
(for  the  same  design  case)  with  the  corresponding  LQG  responses. 

For  the  nonlinear  closed-loop  response,  the  16th-order  compensator  was 
applied  to  Burgers’  equation  and  we  see  (in  Figure  6)  excellent  agreement 
with  the  full  order  closed-loop  trajectory  response.  Hence,  replacing  the 
32nd-order  optimal  LQG  controller  by  a  16th-order  compensator  produces 
a  closed-loop  system  with  minor  performance  degradation. 

We  also  compared  the  performances  of  the  closed-loop  system  of  the  4th- 
order  compensator  with  the  full  order  LQG  responses.  Figures  7  and  8  are 
the  impulse  and  step  responses  of  the  linear  closed-loop  system,  respectively. 
If  one  compares  these  responses  with  the  corresponding  responses  for  the  full 
order  LQG  controller  shown  in  Figures  4  and  5,  then  it  is  clear  that  the  4th- 
order  compensator  performs  almost  as  well  as  the  full  order  LQG  controller. 
Similar  comments  hold  for  the  nonlinear  closed-loop  responses.  For  example, 
the  4th-order  compensator  response  (Figure  9,  solid  line)  compares  well  to 
the  LQG  response  (Figure  9,  Dashed  lines),  especially  after  time  T  —  1.0. 


t 
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Figure  3:  Closed-Loop  TVajectory  for  the  LQG  Case,  N=32 
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Figure  5:  Step  Response  for  the  LQG  &  Fixed-Order 
Cases,  N=32,  Nc=16 


Figure  6:  Closed-Loop  for  the  LQG  &:  Fixed-Order 
Cases,  N=32,  ATC=16 
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Figure  9:  Closed-Loop  for  the  LQ6  &  Fixed-Order 
Cases,  N=32,  Nc—4 


6  Conclusion 

The  purpose  of  this  note  was  to  show  that  finite  dimensional  dynamic  com¬ 
pensators  could  be  used  to  control  a  nonlinear  partial  differential  equation 
without  significant  loss  in  performance.  Although  there  is  considerable  the¬ 
oretical  and  numerical  work  for  bounded  input  and  bounded  output  oper¬ 
ators,  numerical  results  for  the  unbounded  control  and  observation  opera¬ 
tors  are  not  as  fully  developed  as  the  theory.  For  example,  several  authors 
have  considered  questions  of  existence  of  stabilizing  dynamic  compensators 
(even  for  nonlinear  plants  [6])  for  boundary  control  problems.  However, 
approaches,  such  as  the  optimal  projection  method,  that  result  in  a  “com¬ 
putable”  fixed  order  compensator  have  not  been  applied  to  more  general 
boundary  control  problems.  Although  the  numerical  results  presented  here 
show  that  the  optimal  projection  method  can  produce  excellent  designs  for 
problems  with  boundary  control  and  observation,  there  are  a  number  of  the¬ 
oretical  and  numerical  issues  that  need  to  be  resolved  in  order  to  extend  this 
approach  to  practical  problems  of  this  type. 
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